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THE  APPLICATION  OF  FUZZY-SET  THEORY 
TO  THE  BURNTHROUGH  RANGE  EQUATION 


INTRODUCTION 

Electronic  warfare  (EW)  analysis  is  fraught  with  uncertainty,  because  the  tactics  and  systems  of  a 
potential  enemy  are  never  completely  known.  A  continuing  effort  is  made  to  reduce  such  uncertainty 
by  intelligence  collection,  but  there  is  always  a  residue.  Furthermore  the  uncertainty  is  not  statistical. 
If  an  EW  problem  is  treated  by  probability,  then  personal  or  subjective  probability,  that  is,  the  degree  of 
belief  [11,  must  be  used.  Some  authorities  hold  this  to  be  an  improper  use  of  mathematical  probability 
[21.  An  alternative  approach  is  to  use  the  theory  of  fuzzy  sets  [31,  which  is  said  to  be  particularly  suit¬ 
able  for  nonstatistical  uncertainty,  although  agreement  on  this  point  is  incomplete  [4-81.  The  purpose 
of  this  report  is  to  demonstrate  the  application  of  fuzzy-set  theory  to  an  elementary  EW  t  oblem  and 
present  methods  that  may  be  generally  used  with  the  extension  principle  of  fuzzy-set  theory. 

THE  BURNTHROUGH  RANGE  EQUATION 

The  burnthrough  range  equation  for  self-screening  against  a  search  radar  may  be  written  in  the 
form  [91 


*BjEr 

\Kft  08 

4nBREj(Sm)d 

where 

<r  is  the  radar  cross  section  of  the  target  (square  meters), 

Bj  is  the  bandwidth  of  the  jammer, 

Br  is  the  bandwidth  of  the  radar  receiver, 

ER  is  the  effective  radiated  power  of  the  radar,  (ERP)*, 

Ej  is  the  effective  radiated  power  of  the  jammer,  (ERP)  j, 

9W  is  the  radar-antenna  beamwidth, 

9,  is  the  radar-antenna  scan  rate, 

(SNR)rf  is  the  minimum  integrated  signal-to-noise  ratio  necessary  for  detection,  and 
f,  is  the  radar  pulse-rate  frequency. 

Generally,  the  values  of  some  of  the  arguments  in  Eq.  (1)  are  not  precisely  known.  Uncertainty 
in  the  values  of  variables  on  the  right-hand  side  of  Eq.  (1)  is  translated  into  a  corresponding  uncer¬ 
tainty  in  the  value  of  R.  We  apply  fuzzy-set  theory  and  treat  the  uncertain  variables  as  "fuzzy  variables" 
that  may  be  represented  as  a  fuzzy  set,  each  with  a  membership  function.  (The  membership  function 
for  a  fuzzy  set  expresses  on  a  scale  of  0  to  1  the  degree  of  membership  of  each  element  of  some 
universe  of  interest  in  the  set.)  Consider  the  case  where  two  of  the  radar  parameters,  say  ER  and  BR, 
are  uncertain.  Later  we  shall  discuss  the  general  case  where  any  number  of  variables  are  uncertain. 
We  rewrite  Eq.  (1)  as 

RJ-A'2!*-,  (2) 

nR 
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where  K  is  a  known  constant  and  R,  E*,  and  B*  are  the  fuzzy  variables.  Membership  functions  for  E« 
and  BR  are  assumed  given  and  are  denoted  respectively  by  ixE(e)  and  nB(b).  Then  by  Zadeh’s  exten¬ 
sion  principle  UOa).  the  membership  function  for  R  is 

HR(r)  —  suj5  min  (/*£(e),  nB(b)]  (3) 

subject  to  the  condition 

r  -  Kyffb-  (4) 

A  simple  example  of  this  case  with  piecewise  linear  functions  for  nE  and  nB  is  given  in  Appendix  A 
along  with  an  analytical  solution. 

MATHEMATICAL  ASPECTS 

We  first  present  an  analytical  approach  to  the  solution  of  "fuzzy  equations"  using  the  burnthrough 
range  equation  as  an  example. 

Two  Fuzzy  Independent  Variables 


Normalizing  r  in  Eqs.  (3)  and  (4),  we  write 


HR(r)  -  sup  min  (/< E(e),  uB(b) ] 

t,b 

(5) 

subject  to  the  condition 

r2b  —  e. 

(6) 

Substituting  for  e  in  Eq.  (5)  we  obtain 

MB(r)  -  sup  min  \fie(r2b), 

(7) 

We  let  lA.B]  be  the  domain  of  For  each  r  in  both  vr(b)  A  nE(r2b)  and  fiB(b) 

have  the  same  domain,  say  [C,D].  Furthermore,  we  assume  that  both  functions  vr(b)  and  nB(b)  are 
analytic  on  [ C,D\ ■  It  follows  from  the  identity  theorem  [111  that  for  each  r  the  function 
vr(b)  —  fiB(b)  has  at  most  finitely  many  zeros  in  the  interval  \C,D J.  In  other  words,  the  functions 
vr(b)  and  tiB(b)  intersect  at  most  finitely  many  times.  We  let  these  intersections  take  place  at  the 
values  bj,  j  -  1,  2,  ....  rr,  bj  €  [C,£>].  We  now  concentrate  on  a  typical  such  interval,  say  lbj-i,  bj\. 
Without  any  loss  of  generality,  we  assume  that  /ig(b )  lies  below  vr(b)  for  b  €  (by_|,  bj).  Thus,  typi¬ 
cally,  we  have  a  situation  like  the  one  in  Fig.  1. 
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The  points  bJn  »'—  1,  2,  ....  t  in  (b,_ j,  bj)  at  which  dnB/db  —  0  and  simultaneously 
d^fig/db1  <  0  are  the  points  at  which  nB(b)  attains  local  maxima  in  (fy_i,  bj).  We  need  only  examine 

the  points  *,_] ,  bj\,  bj2 . bJt,  bj.  The  largest  of  the  values  /a a(fy_i),  nB(bj,)  with  i  -  1,  2,  ....  t 

and  fiB(bj)  is  equal  to 

sup  mintv,(h),  nB(b)]. 

b 

We  repeat  this  analysis  for  each  of  the  intervals  (fy-i,  bj),  j  —  1.  2,  . . .  ,  n,  using  C  for  b0,  and  for 
the  interval  (b„,  D ),  exchanging  nB  and  v,  when  vr(b)  lies  below  fxB(b).  If  we  denote  the  largest 
such  value  (or  a  largest  value)  in  [bj-\,  A]  by  Mj  and  that  in  [A„,  D]  by  Md,  then  for  the  value  of  r 
used 

M/j(r)  -  max  (Md,  max  Mj).  (8) 

We  then  repeat  the  entire  process  for  enough  values  of  r  to  establish  a  curve  for  n  B. 

Thus  far  we  have  provided  an  analytical  procedure  to  obtain  the  membership  function  for  the 
burnthrough  range  as  a  fuzzy  variable,  provided  that  analytic  membership  functions  for  the  effective 
radiated  power  and  the  bandwidth  of  the  radar  are  given.  The  same  procedure  may  be  used  for  any 
function  of  two  fuzzy  variables  when  one  of  them  can  be  expressed  in  terms  of  the  other  and  the 
dependent  variable,  as  in  Eq.  (6),  and  the  given  membership  functions  are  analytic.  The  last  condition 
is  not  true  for  the  example  of  Appendix  A,  but  local  maxima  are  not  involved  there. 

More  Than  Two  Fuzzy  Independent  Variables 

The  only  theorem  that  seems  to  deal  with  the  more  difficult  case  of  more  than  two  fuzzy  variables 
is  due  to  Baas  and  Kwakemaak  [12, 10b)  and  may  be  stated  as  follows: 

Theorem:  Let  uh  /  -  1,  2 . n,  be  n  piecewise  continuously  differentiable  membership  func¬ 

tions  with  i\.». s  supports.  Let  g  be  a  continuously  differentiable  mapping  of  R"  into  R  (the  real  line). 
At  points  where  the  respective  derivatives  exist,  let  /a ',(*,)  -  dn^xj/dx,  and  £,(x)  - 
dg  (jci,  x2.  ....  x„)/  dx,.  Suppose  that  the  point  x  -  (xl(  x2,  ....  x„)  €  R*  satisfies  the  following: 

•  M'iOc,)  and  g/(x),  /  —  1,2 . n,  exist  and  are  nonzero. 

•  n\(,x\)  “  i*i(x2)  - - //„(x„). 

•  n' ,{xi)l g,(x)  has  the  same  sign  for  each  i  €  [1,  2 . «}. 

Then  x  is  a  strict  relative  maximum  point  of  the  mathematical  programming  problem 

maximize  min  M/(xf) 

/-1.2 . i» 


subject  to  the  condition 

g{Xi,  X2,  ....  X„)  -g(X|,  x2,  ...  ,  x„)  -«(x). 


Suppose  that  values  for  £*,  Bb,  and  fr  in  Eq.  (1)  are  uncertain. 

write 


R  - 


Then  we  define  F  A  if,)0 *  and 


(9) 
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where  k  is  a  constant  and  R,  Ex,  B*,  and  F  are  the  fuzzy  variables.  A  new  membership  function 
H  F(f)  for  F  is  derived  from  a  given  membership  function  for  f,  as  a  fuzzy  variable.  Then  by  the 
extension  principle 


(nR(r)  -  sup  min  1/ iE(e ),  nB(b),  fiF{/)\ 

e.b.j 

subject  to  the  condition 


Normalizing  r  in  Eq.  (9),  we  write 

MR(r)  -  sup  min  \nE(e),  nB(b),  nF(f) J 

t.b.J 

subject  to  the  condition 


r2b  -  ef. 


(10) 


(11) 

(12) 

(13) 


To  apply  the  theorem  of  Baas  and  Kwakernaak,  v/e  must  assume  the  following: 

•  Each  of  the  membership  functions  fiE(e),  fiB(b)  and  nF(f)  is  piecewise 
continuously  differentiable  with  finite  support, 

•  There  exist  e,  b,  and  /such  that  n'E(e),  fi'B(b),  and  n'Fif)  exist  and  are  nonzero; 

•  nE(e)  -  fiB(b)~  nF(/)-, 


A  A  "2 

jji'E(e),  J£n'F(f),  and  ~-~jn'B(b)  have  the  same  sign. 


Under  these  four  conditions  it  follows  from  the  theorem  of  Baas  and  Kwakernaak  that 
HR(r)  is  attained  at  e,  it,  and  j  for  any  given  r.  This  result  assures  us  that  in  principle  under  these 
four  conditions  we  can  find  the  desired  value  of  nR(r).  but  it  does  not  provide  a  procedure.  Further¬ 
more,  the  last  two  conditions  are  rather  strong  and  often  may  not  be  satisfied.  The  theorem  gives  only 
sufficient  conditions  for  a  solution. 


We  now  present  a  method  which  expands  upon  the  previous  one.  This  method  not  only  assures 
the  existence  of  a  solution  under  reasonable  hypotheses  but  also  provides  an  analytical  procedure  for 
finding  fiR(.r). 


We  let  T  - 


so  that  Eq.  (9)  becomes 


R-  WFT. 


(14) 


By  the  procedure  of  the  previous  section  we  can  determine  the  membership  function  /if(r)  for  T;  then 
Eq.  (14)  can  be  viewed  as  an  equation  with  only  two  independent  fuzzy  variables,  F  and  T.  Applying 
the  procedure  once  more,  we  can  determine  the  desired  membership  function  fiR(r)  for  R. 


It  is  important  to  demonstrate  that  we  could  just  as  well  have  chosen  to  make  the  substitution 
T  -  ErF  and  followed  the  same  procedure.  It  is  not  obvious  that  the  substitution  chosen  does  not 
afTed  the  final  result,  that  is,  the  derived  membership  function  nR(r).  That  fiR(r)  is  independent  of 
the  choice  follows  from  a  theorem  dealing  with  the  general  /(-dimensional  case,  which  we  state  and 
prove  in  Appendix  B. 
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Thus  we  have  provided  an  analytical  method  for  obtaining  the  membership  function  for  the 
burnthrough  range,  provided  that  analytic  membership  functions  for  the  effective  radiated  power,  the 
bandwidth,  and  the  pulse  rate  frequency  of  the  radar  are  given.  More  generally,  the  same  procedure 
can  be  used  recursively  when  any  number  of  variables  in  Eq.  (1)  are  fuzzy.  More  generally  still,  the 
method  can  be  used  for  any  fuzzy  equation  where  appropriate  substitutions  can  be  made  and  the  given 
membership  functions  are  analytic. 

If  a  single  "best"  value  of  a  fuzzy  variable  is  desired,  it  may  be  taken  as  the  mean  and  called  the 
(nonstatistical)  expected  value.  Thus  for  burnthrough  range  R,  the  expected  value  is 

£(R)“J*0  rnR(r)  dr/fo  nR(r)  dr.  (15) 

This  definition  is  not  restricted  to  analytic  membership  functions. 


DIGITAL  COMPUTER  SOLUTION 

We  now  discuss  an  approach  which  is  suitable  foT  automatic  computation.  Equation  (4)  can  be 
written  in  the  parametric  form 

e  —  s, 

b  -  ^  (16) 

r 


where  s  is  a  parameter  which  generates  all  points  ( e ,  b)  satisfying  Eq.  (4)  with  fixed  K  and  r.  Substitut¬ 
ing  for  e  and  b  in  Eq.  (5)  gives 


jfJ 

HR(r)  -  sup  min(fz£(s),  /*a(— r $))• 


(17) 


This  is  a  simple  problem  of  maximization  of  a  function  of  a  single  variable,  for  which  standard  numeri¬ 
cal  techniques  are  available.  The  functions  fi£  and  nB  need  not  be  continuous,  much  less  analytic. 


The  method  can  be  applied  generally  to  any  relation  of  the  form 


z  -  fix.  y)  (18) 

for  which  a  suitable  parametric  form 

X  -  Px(z,  s), 

y  -  Py(z,  s)  (19) 

can  be  obtained.  Functions  Px  and  Py  need  not  be  analytic  either,  only  calculable  by  an  algorithm.  The 
membership  function  for  the  dependent  fuzzy  variable  is  then  given  (in  terms  of  membership  functions 
for  the  independent  fuzzy  variables)  by 

#iZ(z)  -  sup  minO/fl/^U  s)l,  n  y[Py(z,  s)l),  (20) 


which  again  only  requires  finding  the  maximum  of  a  function  of  one  variable. 


Figure  2  illustrates  the  method  with  continuous  functions  nx  and  P  y  Loci  of  Eq.  (18)  or,  alter¬ 
natively,  Eqs.  (19),  for  two  values  of  rare  shown  superimposed  upon  contours  of 

Px*r(x,  y)  A  minl/ifOr),  p  y(y)].  (21) 
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Fig.  2  —  Illustration  of  the  computer  solution  for  n2(z)  when  z 
—  f(x.y ),  given  nx(x)  and  ityiy).  Rectangles  represent  con¬ 
tours  of  constant  nXx  y(x.  y)  ~  min  MyO1))  One  set 

of  nested  rectangles  is  obtained  for  convex  nx  and  n  r\  other¬ 
wise,  more  than  one  is  obtained.  Loci  of  /Cx.  y)  are  shown  for 
two  values  of  z.  For  each  value  of  z,  the  corresponding  locus  is 
traversed  by  varying  parameter  s,  and  itx(z)  is  found  as  the 
maximum  of  iix,  r  on  the  locus. 


Since  fix  and  nY  are  continuous  functions,  it  follows  from  Eq.  (21)  that  the  loci  of  constant  nXxY 
(contours)  are  rectangles  (or  degenerate  rectangles,  including  points,  rectangular  areas  or  line  seg¬ 
ments)  with  their  sides  parallel  to  the  x  and  y  axes.  For  each  value  of  z,  the  corresponding  locus  is 
traversed  by  varying  s,  and  /iz(z)  is  found  as  the  maximum  of  nXxY  on  the  locus.  A  FORTRAN  pro¬ 
gram  for  this  procedure  with  the  example  of  Appendix  A  is  presented  in  Appendix  C. 

Equation  (18)  can  be  further  generalized  to  the  case  of  a  function  of  n  fuzzy  variables: 

z  -  /(x,.  x2,  ....  x„).  (22) 

A  parametric  representation  of  /now  has  the  form 
X,  -  P\(z,  5{,  s2.  .... 

X,  -  P„(z,  S|,  s2.  ....  v.|),  (23) 

where  the  S|,  s2 . j„_i  are  independent  parameters.  The  expression  for  fiz  is  then 

#i2(z)  -  sup  min{MilF|(z,  Ji,  s2 . V-i)K  M2lF2(z,  *i.  s2 . •••. 

*i*j . v-i 

M»(F„(z,  5|,  j2,  ... ,  S/f— j)]},  (24) 

where  is  the  membership  function  for  the  fuzzy  variable  X,.  As  the  number  of  parameters 

increases,  numerical  efficiency  becomes  more  important.  In  many  situations  the  membership  functions 
will  be  continuous  and  convex,  allowing  the  use  of  efficient  maximization  algorithms. 

A  nuqor  problem  in  the  application  of  fuzzy-set  theory  is  the  specification  of  membership  func¬ 
tions.  In  the  example  of  Appendix  A,  nE  and  nB  are  assumed  given.  A  more  general  form  of  the  pro¬ 
gram  in  Appendix  C  has  been  written  which  prompts  an  'expert*  user  to  specify  'most  possible,” 
'minimum  possible,”  and  'maximum  possible*  values  of  the  fuzzy  variables.  These  values  are  used  to 
construct  triangular  membership  functions.  A  further  generalization  would  allow  membership  functions 
to  be  specified  for  any  of  the  independent  variables.  Such  a  program  might  be  termed  a  fuzzy  calcula¬ 
tor. 
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DISCUSSION 

The  analytical  and  computer  methods  presented  have  much  wider  application  than  the 
burnthrough  equation  and  EW  analysis.  They  may  be  used  for  any  fuzzy  equation  involving  any 
number  of  variables  when  suitable  substitutions  can  be  made  and,  for  the  analytic  method,  when  the 
given  membership  functions  are  analytic.  The  last  condition  does  not  restrict  the  computer  method, 
which  is  generally  the  most  practical  means  of  solution.  The  analytical  method  is  most  valuable  for 
theoretical  considerations. 

The  expected  advantages  of  the  fuzzy-set  approach  to  analysis  involving  nonstatistical  uncertainty 
are  as  follows: 

•  Membership  functions  (or  agreement  with  suggested  membership  functions)  for  fuzzy 
variables  should  be  easier  to  obtain  from  experts  who  tend  to  be  unwilling  to  commit 
themselves  to  single  numbers. 

•  The  uncertainty  is  made  explicit  and  carried  through  the  analysis  to  appear  as  a  fuzzy 
result.  This  is  less  misleading  than  the  use  of  uncertain,  single  numbers  leading  to  a 
sharp  result. 

•  There  is  a  built-in  sensitivity  analysis  with  simultaneous  variation  of  uncertain  quantities 
and  weighting  of  their  possible  values  (by  means  of  their  membership  functions). 

The  above  advantages  are  cited  against  the  use  of  expected  values  with  a  personal-probability 
(degree-of-belief)  approach.  A  more  proper  comparison  would  be  against  personal  probability  with  dis¬ 
tributions  for,  say,  E*  and  as  random  variables  that  lead  to  a  distribution  for  R  as  a  random  vari¬ 
able.  But  such  an  approach  is  not  currently  used.  The  relative  merits  of  it  and  the  fuzzy-set  method 
need  further  study. 
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Appendix  A 


EXAMPLE  OF  THE  BURNTHROUGH  EQUATION  WITH  TWO  FUZZY 

INDEPENDENT  VARIABLES 


We  assume  the  known  parameters  in  Eq.  (1)  are 

or  —  10s  m2, 

Bj  -  107  Hz, 

Ej  -  106  W, 

(SNR),  -  10, 

0W  *-  ir/180  rad, 

0,  -  ir/5  rad/s,  and 

/,  -  102  Hz. 

We  let  the  membership  (unctions  for  EK  and  BR  as  fuzzy  variables  be  given  by  Figs.  A1  and  A2. 
These  functions  can  te  written 


f  (b  -  106).  b  €  [106,  5-106], 

(Ala) 

-  1  5  b 

4(  5-10* 

J.  ft  €  (510*,  1071, 

(Alb) 

-  1  J  - 
4(  5  10* 

J,  A  €  (107,  I.5-107), 

(Ale) 

and 

(■!<')  “  9  -  >|. 

e  €  (10*,  10’], 

(A2a) 

-  i(*°  - 1?| 

,  *  €  do9,  to10]. 

(A2b) 

It  is  convenient  to  make  the  substitutions  BK  -  Bg -10*, 
b  —  AMO6  with  e  »  *M0*.  Then  Eq.  (2)  becomes 

Er  -  fjJlO®,  and,  correspondingly. 

E  * 

R2-  (10*)J-x, 

(A3) 

which  we  write  as 

R2-  (10*)2T. 

(A4) 

GROSS,  HANRAHAN,  AND  HOOD 


b  (Hi) 


Fig.  Al  —  Membership  function  for  BR 
as  a  fuzzy  variable 


Fig.  A2  —  Membership  function  for  ER 
as  a  fuzzy  variable 
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Our  first  task  is  to  find  the  membership  function  #<r(/)  for  T  with  /  A  e*/b*  (where  A  means  "defined 
as"),  using  the  new  membership  functions: 


-  j(A*  -  1),  b*  €  (1,5], 

(A5a) 

»  1 
“  4 

5-  b •  6  [5,10], 

(A5b) 

3 

"  4 

3  -  y  ,  6*  €  [10,15], 

(A5c) 

and 

1),  e*  €  [1,10], 

(A6a) 

1 

“  9 

r 

10-  ,  *•  €  [10,100], 

1  10 

(A6b) 

According  to  the  extension  principle,  the  membership  function  for  T  is  found  from 

HT(t)  -  sup  min  nl(e*)\, 

b‘.e‘ 

(A7) 

subject  to  /  —  e*/b*,  or  simply 

M  r(/>  •  sup  min  (//*(6*), »»,(£•)], 

(A8) 

where  v,{bm)  A  p l(tb*). 

Figure  A3  is  a  plot  of  na(b')  and  plots  of  fi,(b*)  for  /  —  0.5,  1,  2,  and  3  and  indicates  how  a 
graphical  solution  might  be  obtained.  In  this  simple  case  there  are  no  local  maxima,  and  solutions  are 
found  from  intersections  of  and  i >,(b*).  It  is  also  clear  from  the  figure  which  branches  of  the 

functions  give  a  solution  for  different  regions  of  the  domain  (b*j. 


Thus,  analytically,  we  find 
5/-} 

Mr(t)  -  — i-,  /  €  10.0667,0.775], 

f'+3 


(A9a) 


25/  —  I 
20/  +  9 


.  /  €  10.775,2], 


(A9b) 


100-  / 


90  +  4/ 

0,  otherwise. 


.  t  €  12,100], 


(A9c) 


The  desired  membership  function  for  the  bumthrough  range  is  then  found,  in  view  of  (A4),  from 

HR(r)  -  hh(IOK-Si)  -  nT(t).  (A10) 


This  is  plotted  in  Fig.  A4  as  nR  (r'). 


05  10  15 

r'  (km) 


Fig.  A4  —  Resulting  membership  function 
for  the  burnthrough  range 
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THEOREM  THAT  THE  RESULT  IS  INDEPENDENT 
OF  THE  SUBSTITUTION  USED 


For  a  case  of  n  fuzzy  independent  variables  we  use  the  extension  principle  in  a  more  general 
form*  than  Eq.  (10).  We  let  X  be  a  Cartesian  product  of  universes,  X  —  JIT|  x  •••  x  X„,  and  let 

A  | . A„  be  n  fuzzy  sets  in  A'j,  ....  X„  respectively.  We  let  /be  a  mapping  from  X  lo  a  universe  Y 

such  that  y  -  fix),  y  €  Y,  x  A  (x|,  ...  ,  xn),  x  €  X,  and  x,  €  X,,  i  —  1,  2,  . . .  ,  n.  The  extension 
principle  allows  us  to  induce  from  the  n  fuzzy  sets  A,  with  membership  functions  /ja  Or,)  a  fuzzy  set  B 
on  Y  through  /such  that  the  membership  function  for  B  is 

HBiy)  -  sup  min  l/i.Oc,) . nA  (xs)l 

X  1  " 

y-fix) 

-oifr'o-)-*.  (bi) 


where  /“'O')  is  the  inverse  image  of  y  and  <f>  is  the  empty  set. 

Theorem:  Let  A  be  a  mapping  from  x  Xn  to  a  universe  Z  and  g  be  a  mapping  from 
X\  x  x  X„~2  x  Z  to  the  universe  Y  such  that  for  every  y  €  Y  the  set  of  simultaneous  solutions  (x) 
of  z  -  A(x„_ 1,  x„)  and  y  -  gix{,  ....  x„_2,  z)  is  equal  to  the  solution  set  of  y  -  fix).  Let  C  be  a 
fuzzy  set  induced  on  Z  thorough  h  by  the  fuzzy  sets  A„^  j  and  A„.  Then 

MaO-)-  sup  minl^  (x,) . Ma.  Sxn-3}>  McO)].  (B2) 

*1 . 1 

>-*<*  i . 

Proof:  It  suffices  to  prove  the  theorem  for  n  -  3.  The  general  case  then  follows  by  mathematical 
induction.  Thus  we  want  to  prove  that,  for  all  x  -  (x|,  x2,  x})  in  the  solution  set  of  y  —  fix), 

sup  mink^.UO.At^jUjJ.M^^xj))  - 
Z-/U) 

sup  minkx,k|),  sup  mink,4,Or2),M,r,(*j)l}-  (B3) 

*|.l  1  Xj.JTj  1  3 

y-*(x,.z)  j-a(xj,x3) 


We  let 

S/O')  -  |x:/(x)  -  y.  mink.x,(*i)-M/42(x2),M/«jOrj)]  -  ii^ix,)).  (B4) 

3 

where  i  —  1,  2,  3.  It  is  clear  that  S/O')  is  the  solution  set  of  y  «  fix).  It  suffices  then  to  show 

/-i 

that  Eq.  (B3)  is  valid  for  all  S/O')- 


*D.  Dubois  and  H.  Pride,  Fuzzy  Sets  and  Systems:  Theory  and  Applications,  Academic  Press,  New  York,  1980,  pp.  36-37. 
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We  first  evaluate  the  right-hand  side  of  Eq.  (B3)  for  Si  O’).  From  definition  (B4)  we  have  for  all 

x  es,0> 


Ma,(* i>  <  minfo^x^./i^xj)! 

<  sup  minl/r^  (x2),ma ,Oj)1- 

Jtj.Xj 

Thus  the  right-hand  side  of  Eq.  (B5)  becomes  for  x  €  St(y) 

sup  mink 4, Or,),  sup  mink^Or^./^/xj)]}  -  sup  fiA  (x,). 

X,.Z  1  Xj.X3  1  3  X, 

2—h(X2,Xy) 

But  for  x  €  Siiy)  the  left-hand  side  Eq.  (B5)  can  be  written,  in  view  of  Eq.  (B4),  as 

sup  minfe  ,  (xi).>i .  (jr2).M.4,(*a))  -  sup/^.Or,). 

X  1  *  *  Xj  ' 

Thus  Eq.  (B5)  is  valid  for  x  €  S]0)- 

For  S20)  we  have  from  Eq.  (B4) 

Ma2(x2)  -  mink^Or^.M^Orj)), 
so  that  in  the  right-hand  side  Eq.  (B3)  for  x  €  S20) 

sup  minkx,^"^,^)]  “  supm,4,(x2). 


(B5) 


(B6) 


(B7) 


(B8) 


(B9) 


Xj.Xj  -  '  <2 

Also,  since  /i/4j(x2)  <  Ha^x\)  f°r  S20)»  it  follows  that  the  right-hand  side  of  Eq.  (B3)  becomes  for 
x  €  S20> 


sup  mink^.Oi).  sup  Ma,(*2»  “  sup  /a a ,0r2). 


x,.l 


*2 

*«-A(x  j.Xj) 


For  S20)  the  left-hand  side  of  Eq.  (B3)  becomes  with  (B4) 

sup  “  SUP  Ha^xJ- 

Thus  Eq.  (B5)  is  valid  for  x  €  S20)  as  well. 


(BIO) 


(BID 


Similarly  we  can  show  that  Eq.  (B5)  is  valid  for  x  6  S3(y)  and,  therefore,  for  all  x  in  the  solution 
set.  The  hypotheses  of  the  theorem  are  sufficient;  we  can  show  that  they  are  not  necessary.  Further, 
the  hypotheses  are  satisfied  with 

/Or, . x2)  -  *(x, . x,_2.Mx„_,,x,,)J.  (B12) 

The  problem  of  when  /  has  the  desired  property  is  interesting  and  will  be  discussed  in  a  subsequent 
report. 
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Appendix  C 

COMPUTER  SOLUTION  OF  THE  EXAMPLE  IN  APPENDIX  A 

The  FORTRAN  program  listed  in  Table  Cl  calculates  the  membership  function  and  average 
(mean)  value  of  the  bumthrough  range  for  the  example  in  Appendix  A.  To  show  the  generality  of  the 
procedure,  the  parametric  expressions  for  e  and  b  and  the  functions  nE  and  ixg  are  implemented  as 
FUNCTION  subprograms.  The  subroutine  LIMITS  restricts  the  search  to  regions  of  the  ( e ,  b )  plane 
where  n  E  and  /a  B  are  nonzero.  Thus  the  program  is  applicable  to  any  case  where  the  effective  radiated 
power  and  the  bandwidth  are  the  only  uncertain  radar  parameters.  Results  are  shown  in  Table  C2  for 
the  membership  functions  and  known  parameter  values  of  Appendix  A. 


Table  Cl  —  Computer  program  for  the  burn-through  range  when  the 
radar’s  effective  radiated  power  and  bandwidth  are  uncertain 


PROGRAM  FUZZ 

c 

(^* ******************************* ************************************** 

C  A  PROGRAM  FOR  CALCULATING  THE  BURNTHROUGH  RANGE  FOR  A  HOSTILE 

C  SEARCH  RADAR  AGAINST  A  SELF-SCREENING  JAMMER. 

C 

0********************************************************************** 

C 

c* ************************. **.***•****.***•*******.**.******. ********** 

c 

C  THE  BURNTHROUGH  RANGE,  R,  DEPENDS  ON  THE  ERP,  E,  AND  BANDWIDTH, 

C  B,  OF  THE  RADAR  THROUGH:  R**2=A*B/B,  WHERE  A  DEPENDS  ON  KNOWN 

C  PROPERTIES  OF  THE  JAMMER  AND  THE  RADAR. 

C  THIS  PROGRAM  CALCULATES  THE  MEMBERSHIP  FUNCTION  FOR  THE  FUZZY 

C  VARIABLE  R  WHICH  DEPENDS  ON  THE  FUZZY  VARIABLES  E  AND  B. 

C  FOR  GIVEN  R  WE  HAVE  E«PE(R,S),  B*PB(R,S),  WHERE  S 

C  IS  A  PARAMETER  WHICH  GENERATES  ALL  E,B  PAIRS  WHICH 

C  YIELD  R.  E  AND  B  HAVE  MEMBERSHIP  FUNCTIONS  FMUE(E)  AND  FMUB(B). 

C 

c* ******************** ************************************************* 

C 

PARAMETER  PI=3. 14159 
COMMON  SMIN , SMAX , RMIN , RMAX , A 
NR=100 
NS-1000 
C 

TYPE  * , 1  RADAR  CROSS  SECTION  (SQ.  M)  ?' 

ACCEPT  *,SIG 

TYPE  JAMMER  BANDWIDTH  (MHZ)  ?' 

ACCEPT  *,BJ 

TYPE  *, 'JAMMER  ERP  (WATT)  ?' 

ACCEPT  *, EJ 

TYPE  RADAR  BEAMWIDTH  (DEG)  ?' 

ACCEPT  *,W 
W»W*P 1/180. 

TYPE  * , ' SCAN  RATE  (HZ)  ?' 

ACCEPT  *,SR 
SR»2.*PI*SR 

TYPE  * , ' MIN .  SNR  FOR  DETECTION  ?’ 

ACCEPT  *,SNR 

TYPE  • , ' PRF  (PPS)  ?' 

ACCEPT  * , PRF 
C 


Table  continues 
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Table  Cl  (continued) 


A=  (SIG*BJ/(4  ,‘PI‘EJ‘SNR)  )  *  (W*PRF/SR)  “0 . 8 
A*A‘1E8 
C 

CALL  LIMITS (0 . ,0) 

SUM-0  , 

WSUM-0  • 

AVG=0. 

C 

DO  500  1*1, NR 

R=RMIN+  (RMAX-RMIN )  ‘FLOAT  (1-1 )  /FLOAT  (NR-1 ) 

CALL  LIMITS (R,l) 

SUP*0. 

DO  400  J*1,NS 

S=SMIN+  (SMAX-SMIN)  ‘FLOAT  (J-l )  /FLOAT  (NS-1 ) 
SUP=AMAX1 (SUP, AMIN  1 (FMUEfPE (R,S) ) ,FMUB(PB (R,S) ) ) ) 
400  CONTINUE 

R=R/1000. 

TYPE  *,'R  (KM)  =  ',R,‘MU(R)  =',SUP 
C 

C  CALCULATE  EXPECTED  VALUE 
C 

SUM=SUM+SUP*R 

WSUM-WSUM+SUP 

C 

500  CONTINUE 

C 

IF  (WSUM.GT.lE-6)  AVG>SUM/WSUM 

TYPE  ‘.‘AVERAGE  BURNTHROUGH  RANGE  (KM)  »‘,AVG 


c$$ $$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
c 

SUBROUTINE  LIMITS (R, I) 

C 

c  * 

c  THIS  SUBROUTINE  DETERMINES  SEARCH  LIMITS  ON  R  AND  S.  * 

C  IT  REQUIRES  UPPER  AND  LOWER  LIMITS  ON  E  AND  B  * 

C  OUTSIDE  OF  WHICH  FMUE(E)  AND  FMUB(B)  ARE  ZERO,  * 

C  AND  EXPRESSIONS  FOR  R  AND  S  IN  TERMS  OF  E  AND  B.  * 

C  * 

c 

COMMON  SMIN , SMAX, RMIN , RMAX , A 
C 

DATA  EM  IN/1 .  / ,  EMAX/1 00 .  / ,  BM  IN  /I .  / ,  BMAX/1 5 .  / 

C 

IF  (I.GT.0)  GOTO  100 
RMIN-SQRT (A*EMIN/BMAX ) 

RMAX-SQRT (A*EMAX/BMIN ) 

RETURN 

C 

100  CONTINUE 

SMIN-AMIN1  (EMIN,BMIN*R“2/A) 

SMAX-AMAX1 ( EM AX , BMAX  *R  “  2 /A ) 

RETURN 


Table  continues 


16 


NRL  REPORT  8619 


FUNCTION  PE(R,S) 


THIS  FUNCTION,  TOGETHER  WITH  PB(R,S>,  DEFINES  THE  DEPENDENCE 
OF  R  ON  VARIABLES  E  AND  B.  THE  RELATION  IS  R**2«A*E/B. 


COMMON  SMIN , SMAX , RMIN , RMAX , A 

C 

PE-S 

RETURN 

END 

C 

C$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

C 

FUNCTION  PB (R,S) 

C 


C  THIS  FUNCTION,  TOGETHER  WITH  PE(R,S),  DEFINES  THE  DEPENDENCE  * 

C  OF  R  ON  VARIABLES  E  AND  B.  THE  RELATION  IS  R*»2-A»E/B.  * 

C  * 

>*»*»**»*«*.»**»»»*» 

C 

COMMON  SMIN, SMAX, RMIN, RMAX, A 

C 

PB-A*S 

IF  (R.CT.1E-3)  PB-PB/R**2 

RETURN 

END 

C 

0$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

C 

FUNCTION  PMUE(E) 

C 


THIS  REPRESENTS  THE  MEMBERSHIP  FUNCTION  FOR  FUZZY  VARIABLE  E. 


C  SIMPLE  EXAMPLE 
FMUE-0. 

IF  (E.GE.1..AND.E.LE.10.)  FMUE- (E-I .1/9. 

IP  (E.GT. 10. . AND. E. LE. lt«. )  FWUE- <1 0 . -E/1 0 . ) /9 . 
RETURN 
END 
C 
C 

FUNCTION  FMUB(B) 


THIS  REPRESENTS  THE  MEMBERSHIP  FUNCTION  FOR  FUZZY  VARIABLE  B. 


SIMPLE  EXAMPLE 
FMUB “0  - 

IF  (B.GE. 1 . . AND. B. LE . 5. )  FMUB- (B-l . ) /4 . 

IF  (B.GT.5. .AND.B.LE.10.)  FMUB- (5  . -B/5 . ) /4  . 

IF  (B.GT.10. .AND.B.LE.15.)  FMUB- (3 . -B/5 . ) *3 . /4 . 

RETURN 

END 
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Table  C2  —  Computer  output  for  the  program  in  Table  Cl  with 
the  membership  functions  and  known  parameter  values  of 
Appendix  A 


*  RUN  FUZZ 

RADAR  CROSS  SECTION  (S 
1.0ES 

JAMMER  BANDUIDTH  (MHZ) 
10.0 

JAMMER  ERP  ( IJATT )  ? 

1 .0E6 

RADAR  BEAM  UIDTH  (DEG) 
1.0 

SCAN  RATE  (HZ)  ? 

0.1 

MIN.  SNR  FOR  DETECTION 
10.0 

PRF  (PPS)  ? 

100.0 

R.(K'M)  *  0.3466001 

0.  M)  ? 

? 

? 

? 

MU(R ) 

* 

0 . 0000000E+00 

R  (KM) 

s 

0.4786926 

MU  ( R ) 

= 

9 . 185S764E-0H 

R  (KM) 

s 

0.6107852 

MU(  R ) 

= 

0.2004713 

R  (KM) 

* 

0.7428777 

MU(R ) 

5 

0.3220916 

R  (KM) 

s 

0.8749703 

MU(R ) 

■ 

0.4459961 

R  (KM) 

9 

1.007063 

MU(R ) 

X 

0.5818230 

R  (KM) 

= 

1.139155 

MU(R ) 

■ 

0.7087278 

R  (KM) 

9 

1.271248 

MU(P) 

s 

0.7931786 

R  (KM) 

9 

1.403341 

MU(R ) 

9 

0.8506137 

R  (KM) 

• 

1.535433 

MU(R ) 

9 

0.9012267 

R  (KM) 

m 

1.667526 

MU(R ) 

9 

0.9414500 

R  (KM) 

s 

1.799618 

MU(R) 

S 

0.9768122 

R  (KM) 

M 

1.931711 

MU(R ) 

9 

0.9953954 

R  (KM) 

• 

2.063803 

MU(R) 

9 

0.9810811 

R  (KM) 

• 

2. 195896 

MU(R ) 

9 

0.9656657 

R  (KM) 

9 

2.327989 

MU(R ) 

9 

0.9502503 

R  (KM  > 

* 

2.460081 

MU(R) 

9 

0.9337337 

R  (KM) 

• 

2.592174 

MU(R ) 

9 

0.9172173 

R  (KM) 

V 

2.724266 

MU(R ) 

9 

0.8995996 

R  (KM) 

• 

2.856359 

MU(R ) 

9 

0.8830831 

R  (KM) 

9 

2.988451 

MU(R ) 

9 

0.8651874 

R  (KM) 

« 

3.120544 

MU(R ) 

9 

0.8467467 

R  (KM) 

m 

3.252636 

MU(R ) 

9 

0.8291292 

R  (KM) 

9 

3.384729 

MU(R ) 

9 

0.8104104 

R  (KM) 

9 

3.516821 

MU(R ) 

9 

0.7926978 

R  (KM) 

9 

3.648914 

MU(R) 

9 

0.7738354 

R  (KM) 

9 

3.781007 

HU(R ) 

9 

0.7561366 

R  (KM) 

9 

3.913099 

MUCR) 

9 

0.737109? 

R  (KM) 

9 

4.045192 

MU(R) 

9 

0.7195187 

R  (KM) 

9 

4.177284 

MU(R ) 

9 

0.7004953 

R  (KM) 

9 

4.309377 

MU(R) 

9 

0.6831930 

R  (KM) 

9 

4.441469 

MU(R ) 

9 

0.6645205 

R  (KM) 

9 

4.573562 

MU(R ) 

9 

0.6475086 

R  (KM) 

9 

4.705654 

MU(R ) 

9 

0.6290241 

R  (KM) 

9 

4.837747 

MU(R  i 

= 

0.6128144 

R  (KM) 

9 

4.969840 

MU(R) 

■ 

0.5948097 

R  (KM) 

9 

5.101932 

MU(R ) 

■ 

0.5781121 

R  (KM) 

9 

5.234025 

MU(R ) 

9 

0.56213 02 

R  (KM) 

9 

5.36611? 

MUCR) 

9 

0.5451371 

R  (KM) 

9 

5.498210 

MU(R) 

9 

0.5235189 

Table  continues 
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Table  C2  (continued) 


R 

(KR) 

X 

5.630302 

RU(R ) 

_ 

0.5123141 

R 

OCR) 

X 

5.762395 

MU(R> 

a 

0.4966409 

R 

(KR ) 

a 

5.894488 

RU(R) 

a 

0.4813002 

R 

(KR) 

X 

6.026580 

RU(R ) 

= 

0.4665687 

R 

OCR) 

s 

6.158673 

RU(R ) 

s 

0.4524058 

R 

(KR) 

X 

6.290765 

RU(R) 

a 

0.4387506 

R 

(KR) 

X 

6.422858 

RU(R  ) 

s 

0.4233998 

R 

(KR) 

% 

6.5S4951 

RU(R ) 

= 

0.4097492 

R 

(KR) 

M 

6.687043 

RU(  R ) 

a 

0.3964981 

R 

(KR) 

X 

6.819135 

RU(R) 

a 

0.3828891 

R 

(KR) 

9 

6.951228 

RU(R ) 

a 

0.3696639 

R 

(KR) 

m 

7.083320 

RU(R ) 

c 

0.3581413 

R 

(KR) 

a 

7.215413 

RUCR) 

a 

0.3449887 

R 

(KR) 

9 

7.347506 

RU(R ) 

a 

0.3319885 

R 

(KR) 

n 

7.479598 

RU(R) 

« 

0.3199065 

R 

(KR) 

m 

7.611691 

RU(R) 

a 

0.3084331 

R 

(KR) 

m 

7.743783 

RU(R) 

a 

0.2969701 

R 

(KR) 

X 

7.875876 

R(J(R ) 

a 

0.2860853 

R 

(KR) 

9 

8.007968 

RU(R ) 

a 

0.2762794 

R 

(KR) 

a 

8.140061 

RU(R) 

s 

0.2663836 

R 

(KR ) 

9 

8.272153 

RU(R ) 

a 

0.2549562 

R 

(KR) 

9 

8.404246 

RU(R ) 

a 

0.2448O68 

R 

(KR) 

W 

8.536339 

RU(R ) 

a 

0.2358705 

R 

(KR) 

9 

8.668430 

RU(R ) 

a 

0.2248322 

R 

(KR) 

9 

8.800524 

RU(R) 

a 

0.2168122 

R 

(KR) 

9 

8.932616 

RU(R) 

a 

0.2074694 

R 

(KR) 

9 

9.064708 

HU(R ) 

a 

0.1977779 

R 

(KR) 

X 

9.196801 

RU(R ) 

s 

0.1901431 

R 

(KR) 

9 

9.328894 

RU(R) 

a 

0. 1825 137 

R 

(KR) 

r 

9 . 460986 

RU(R ) 

a 

0.1730798 

R 

(KR) 

9 

9.593080 

RU(R ) 

a 

0.1639993 

R 

(KR) 

a 

9.725171 

RU(R ) 

s 

O.1559047 

R 

(KR) 

s 

9.857264 

RU(R ) 

a 

0.1482931 

R 

(KR) 

a 

9.989355 

RU(R) 

s 

0.14O6857 

R 

(KR) 

9 

10.12145 

RU(R) 

a 

0.1337029 

R 

(KR) 

a 

10.25354 

Rll(R ) 

a 

0.1277320 

R 

(KR) 

• 

10.38563 

RU(R ) 

a 

0.1216348 

R 

(KR) 

a 

10.51773 

RU(R) 

a 

0.1140414 

R 

(KR) 

a 

10.64982 

RU(R) 

a 

0.1064510 

R 

(KR) 

» 

10.78191 

RU(R) 

a 

0.1000918 

R 

(KR) 

a 

10.91400 

RU(R ) 

a 

9 . 5029056E-02 

R 

(KR) 

a 

11.04610 

RU(R ) 

a 

8 . 74472E6E-02 

R 

(KR) 

u 

11.17819 

RU(R ) 

a 

8.2949318E-02 

R 

(KR) 

9 

11.31028 

RU(R ) 

a 

7.60412S1E-02 

R 

(KR) 

9 

11.44238 

RUCR) 

a 

7.0592560E-02 

R 

(KR) 

9 

11.57447 

RUCR) 

a 

6.4644873E-02 

R 

(KR) 

9 

11.70656 

RUCR ) 

a 

6.0495801E-02 

R 

(KR) 

m 

11.83865 

RUCR) 

a 

5.3257197E-02 

R 

(KR) 

9 

11.97075 

RU(R) 

a 

4.9441695E-02 

R 

(KR) 

9 

18.10284 

RUCR) 

a 

4.3067295E-02 

R 

(KR) 

9 

18.23493 

RU(R) 

a 

3.8065791E-02 

R 

(KR) 

9 

18.36702 

RU(R ) 

a 

3 . 42S5862E-02 

R 

(KR) 

9 

12.49911 

RU(R ) 

a 

3.0447721E-02 

R 

(KR) 

- 

12.63121 

RUCR) 

a 

2.2890359E-02 

R 

(KR) 

9 

12.76330 

RUCR) 

a 

1.9085437E-02 

R 

(KR) 

a 

12.89539 

RUCR) 

a 

1.5282005E-02 

R 

(KR) 

a 

13.08748 

RUCR) 

a 

1 . 1480063E-02 

R 

(KR) 

a 

13.15958 

RUCR) 

a 

7.6795518E-03 

R 

(KR) 

a 

13.29167 

RUCR) 

a 

3.8804412E-03 

R 

(KR) 

a 

13.42376 

RUCR) 

a 

0.0000000E+00 

AVERAGE 

t 

BURNTHROUGH 

RANGE  (KR) 

a 

4.430007 

END 


DATE 

FILMED 


DTIC 


